NMR measurements of hyperpolarized 3 He gas diffusion 
in high porosity silica aerogels 
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Hyperpolarized 3 He is used to nondestructively probe by NMR the structure of custom-made 
and commercial silica aerogels (97% and 98.5% porous). Large spin-echo signals are obtained at 
room temperature and very low magnetic field (2 mT) even with small amounts of gas. Attenuation 
induced by applied field gradients results from the combined effects of gas diffusion and confinement 
by the porous medium on atomic motion. Nitrogen is used as a buffer gas to reach equivalent 
3 He pressures ranging from 5 mbars to 3.5 bars. The observed pressure dependence suggests a 
non-uniform structure of the aerogels on length scales up to tens of micrometers. A description by 
broad phenomenological distributions of mean free paths is proposed, and quantitatively discussed 
by comparison to numerical calculations. The investigated aerogel samples exhibit different effective 
diffusion characteristics despite comparable nominal porosities. 

Keywords: Silica aerogel; spin diffusion; hyperpolarized helium; low field NMR. 



Article published in the Journal of Chemical Physics. 
Reference : J. Chem. Phys. 123, 064506 (2005) 



I. INTRODUCTION 

Silica aerogels are highly porous solids obtained by sol- 
gel process [1]. Their very open structure consists of a 
delicate, highly ramified backbone network of randomly 
interconnected thin silica strands, formed by aggrega- 
tion of small silica clusters [2, 3]. Mechanically robust, 
they exhibit very low densities, thermal conductivities, 
dielectric constants and sound velocities, as well as poor 
chemical reactivity and large specific areas. Despite high 
synthesizing costs and difficulties, they have been com- 
mercially produced for civil and military use and they 
have found applications in thermal insulation, manufac- 
ture of materials and coatings, catalysis, gas storage and 
filtering, Cerenkov radiation detection, etc. In science, 
their unique properties make them attractive for cutting- 
edge studies in optics, acoustics, chemical physics, low- 
temperature physics, nuclear physics, space investiga- 
tions, laser tests, microelectronics, electrical engineering, 
life sciences and bioengineering [1, 4, 5]. The density and 
microstructure of the aerogels can be precisely controlled 
by adjusting the growth chemistry [4, 6-8]. This makes 
them particularly suited for a variety of investigations 
on bulk, film or layer systems in weakly confined random 
geometries. 

Considerable effort has been devoted to the study of 
phase transitions in aerogel [9-12], in particular with he- 
lium at low temperature to probe statistical and quantum 
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effects [10, 12, 13]. The silica network provides a struc- 
tural disorder background to the fluid where it is imbed- 
ded, but occupies only a very small fraction of the sample 
volume. Recently, extremely high (up to 99.5%) porosity 
aerogels have been used to study impurity scattering ef- 
fects on unconventional pairing states in superfluid 3 He 
[14, 15] and changes of critical exponents in superfluid 
4 He [16]. Strong modifications of the phase diagram oc- 
cur when the correlation length of the order parameter 
in the fluid becomes comparable to that of aerogel. The 
microstructure of the samples used in the experiments 
must hence be characterized at the appropriate length 
scales, that range from nanometers (15—80 nm in super- 
fluid 3 He [17], 3-100 nm in superfluid 4 He [16, 18]) to 
tens of microns (mean free path, mfp, of unpaired 3 He 
quasiparticles at low pressure and temperature [12]). 

In contrast with Vycor glass or porous gold, for in- 
stance, there is no single or dominant length scale to 
characterize the distribution of pores and silica in aero- 
gel. Aerogel is known to be fractal at very small scales, 
a feature revealed both by small angle scattering (SAS) 
measurements with x rays or neutrons, and computer 
simulations of the growth process [2, 19, 20]. Self-similar 
structures range over correlation length scales from a few 
nanometers to hundreds of nanometers (depending on 
catalysis conditions and kinetics of the gelation process), 
above which the medium appears to be fairly homoge- 
neous. Actual homogeneity at microscopic level, and rel- 
evance of a description solely based on the fractal dimen- 
sion and one correlation length parameter, are key issues 
for the interpretation of all physical measurements [12], 
acutely debated in the case of superfluid 3 He and 4 He 
in aerogel [15, 21]. Discrepancies between experimental 
results are common and, most generally, tentatively at- 
tributed to unknown sample-to-sample microscopic vari- 
ations (due to growth conditions, or insufficient charac- 
terization [12]). In addition, meso- or macroscopic in- 
homogeneities may develop during sample preparation, 
mounting, or use at low temperature (due to mechani- 
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cal or thermal stress [22]). Unfortunately, conventional 
techniques such as mercury intrusion porometry, ther- 
moporometry, or nitrogen adsorption/desorption, are not 
suitable to characterize the aerogel pore structure accu- 
rately [23]. 

Nuclear magnetic resonance (NMR) is an established 
noninvasive technique to characterize porous media on 
a microscopic scale [24, 25]. Gas probes are of particu- 
lar interest, because they can easily penetrate the porous 
matrix and exhibit fast diffusion rates that are conve- 
niently varied over a wide range of values by pressure 
control. Compared to other inert species, 3 He may be 
considered as an ideal candidate probe for diffusion mea- 
surements: this spin-1/2 atom has long intrinsic relax- 
ation times, and a resonance frequency fairly insensitive 
to most changes of its local environment [26, 27]. Suit- 
able sensitivity is usually achieved using both high pres- 
sure and high magnetic field [27, 28] to compensate for 
the low density of the gas. With laser optical pump- 
ing techniques, high nuclear polarizations (hyperpolar- 
izations) can be achieved independently of the magnetic 
field strength, and very high magnetization densities are 
obtained even at low pressures. Submillimetric imaging 
and diffusion measurements of hyperpolarized gases in 
free or confined spaces have been demonstrated both at 
high and low magnetic field [29-32], and used for various 
purposes including in vivo characterization of the lung 
alveolar airspaces for medical applications [33-36] . In the 
present work, we use hyperpolarized 3 Hc NMR to non- 
destructively probe the void spaces in high porosity silica 
aerogels identical to those used in most low-temperature 
experiments. With probe gas mfp's ranging from a frac- 
tion of micrometer to a few tens of micrometers when 
pressure is varied, this technique yields a unique type of 
information on the structure of these materials. 

We had previously performed preliminary measure- 
ments on aerogel samples to test the feasibility of room 
temperature NMR investigations [37-39]. They sug- 
gested a non-uniform distribution of density of silica 
strands in this system, in contrast, for instance, with the 
single quasiparticle mfp value obtained by spin diffusion 
measurements in 3 He- 4 He dilute solutions [40]. We have 
now undertaken a more extended investigation to con- 
firm this result, using an improved experimental setup. 
In this article, we report on a systematic study over a 
broader range of gas pressures that has been carried out 
to probe several samples of custom-made and commer- 
cially produced aerogels. Our former conclusions [38, 39] 
are confirmed. We show that all aerogels exhibit the same 
qualitative behavior, characterized by a pressure depen- 
dence of the gas diffusion coefficient that excludes a de- 
scription of the aerogel by a single characteristic length. 
Significant differences are observed between the various 
samples, that do not correlate with the specified aerogel 
porosities and may rather be associated to differences in 
growth conditions. We also propose a phenomcnologi- 
cal description of gas diffusion in aerogel based on broad 
power-law distributions of mfp's. Their main features are 



quantitatively derived from the NMR data analysis, and 
discussed on the basis of computer calculations. 



II. EXPERIMENT 
A. Aerogel samples 

Six samples of aerogel with 97% and 98.5% nominal 
porosities have been studied. One is custom made, and 
all others come from specialized manufacturers. 

Sample M is a 98.5% open custom-made aerogel 
prepared by hydrolysis and polymerization of tctram- 
ethoxysilane, dilution in acetonitrile, gelation after addi- 
tion of ammonia, and supercritical drying in an excess of 
solvent to preserve the delicate structure [41]. Because 
of the base catalyzing, the gelation process is believed 
to be diffusion-limited aggregation of small silica clus- 
ters (1—3 nm in diameter). This type of aerogel, used 
in many low-temperature experiments, has been fairly 
characterized at small length scales using experimental 
(scanning and transmission electron micrography, SEM 
and TEM; SAS of x rays and neutrons) and numeri- 
cal (diffusion-limited cluster aggregation, DLCA, growth 
models) techniques [20] . The nanostructure of computer- 
generated aerogels, nicely reproducing that of real sam- 
ples probed by SEM/TEM and SAS, can be most con- 
veniently visualized (see, e.g., [12, 20] for 2D or 3D im- 
ages) and statistically analyzed. On custom-made 98.5% 
samples similar to ours, SAS of x rays shows a fractal- 
like distribution of silica extending up to a length scale of 
65—100 nm [13, 16], from which a purely geometrical mfp 
on the order of 270 nm can be calculated [42] . In a com- 
parable DLCA model (98.2% porosity), the very dilute 
aerogel structure is found to be organized so that half 
of the widely open volume lies farther than 10 nm from 
some silica, while none of it is farther than 35 nm [42]. 
Recently, path-length distributions have been analyzed 
in simulated aerogels using randomly sampling straight- 
line trajectories [43], and finite-size hard spheres dis- 
placements to mimic experimental conditions [44]. The 
estimated 3 He quasiparticle mfp significantly exceeds in 
aerogel that of a randomly distributed medium of iden- 
tical volume density (180 and 120 nm, respectively, for a 
98.5% open sample [44]). 

Aerogel M has been grown in a cylindrical glass con- 
tainer (inner diameter 12 mm, height 37 mm) partly filled 
with a set of three glass tubes (diameter 6x4 mm, height 
30 mm, and axes parallel to that of the container). It 
has been produced as a transparent, crack-free sample 
attached to the glass walls. An epoxy end cap for the 
aerogel container is glued to allow leak-tight connection 
to the 3 He polarization system, with negligible dead vol- 
ume between the free aerogel surface and the epoxy cap 
(Stycast 1266 A/B, Emerson and Cuming, Billerica, MA, 
USA). 

Sample A is a 97% porous aerogel with nominal specific 
area 465 m 2 /g manufactured by Airglass, Staffanstrop, 
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Sweden. The cylindrical sample (diameter 12.5 mm, 
height 10.8 mm) is enclosed in a matching PMMA [45] 
container. The container and end cap are machined so 
that the aerogel gently fits in, avoiding mechanical stress 
(^0.05 mm gap size). 

Samples Ji— J4 are four pieces of 97% porous aero- 
gel with nominal specific area 510 m 2 /g manufactured 
by Matsushita Electric Works, Osaka, Japan. Ji and 
J2 originally had identical parallelepipedic shapes (sec- 
tions 15x15.5 mm 2 , thicknesses 9.7 mm). Ji, acci- 
dentally crushed before use, is reduced in thickness by 
1.7—2.7 mm (from end to end along its longest dimen- 
sion, with a noticeable transverse fracture line). J3 and 
J 4 have cylindrical shapes (diameters 9.1 and 7.7 mm, 
lengths 10.7 and 10.5 mm, respectively). Their original 
rough ends have been cut flat using a silicon carbide disk 
driven at very high speed (65 000 rpm). Viewed through 
these new flat ends, the aerogels are perfectly transpar- 
ent and display no macroscopic flaw. All four samples 
are enclosed in tightly fitting PMMA containers (apart 
from a 0—1 mm gap left open for the trapezoidal sample 
Ji). 

An empty PMMA container, identical to that of aero- 
gel A (diameter 12.5 mm, height 10.8 mm), is used as a 
reference cell for free diffusion measurements. 

Polypropylene Luer fittings are glued to the PMMA 
containers to allow connection to the 3 He polarization 
system, through a 10-cm-long Tygon tube (Masterflex, 
Cole-Parmer Instrument Co., Vernon Hills IL, USA) with 
1.6-mm inner diameter. A 0.7xl.l-mm Teflon tube is in- 
serted into the Tygon filling tube to reduce the dead vol- 
ume. This also increases the gas flow impedance, hence 
limits the pressure change rates inside the aerogel sample. 

There has been no attempt to clean the aerogels. The 
surface may have retained traces of preparation chemi- 
cals, as well as O2 or moisture adsorbed during storage. 
During experiments, the samples are under vacuum, or 
in contact with high purity 3 He or N2 gas. A totally 
negligible amount of He atoms may dissolve in the silica 
[27]. The aerogel behavior is controlled during and after 
use by visual inspection. Sample M remains clearly well 
attached to the glass walls. Several cracks are present 
in its bulk part, between the free surface and the top 
of the inner tubes that presumably provide mechanical 
support to the aerogel. They developed during prior use, 
due to rough filling and emptying procedures [37], but 
remained unaffected both by the breaking of one inner 
tube and the present scries of measurements. All off- 
the-shelf samples (A, Ji— J4) are used for the first time. 
Through the PMMA side walls, Ji and J2 can be ob- 
served to expand and shrink slightly as pressure varies 
from to 1 atm. Very few cracks developed inside these 
commercial aerogels, allowed to move freely inside the 
containers by an irreversible decrease of all their dimen- 
sions (by about 5%— 10%) observed after a few pressure 
cycles. 



B. NMR setup 

The pulsed NMR experiments are performed at low 
frequency (62.4 kHz) using a homemade system, in the 
homogeneous static field used for optical pumping (OP) 
described in Section II C. We use free induction de- 
cay (FID) signals to monitor the 3 He magnetization, 
and spin-echoes generated by a standard Carr-Purcell- 
Meiboom-Gill (CPMG) sequence with an applied gradi- 
ent to measure diffusion [46]. 

All NMR and gradient coils depicted in Fig. 1 are cop- 
per wire windings held by PMMA stands. They are sur- 
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FIG. 1: Sketch of the NMR coils. G: set of four gradient coils 
(45x42 (45x35)-mm size, 4-mm horizontal width, 38 (78)-mm 
average spacing, and 9 (25) turns each, for inner (outer) coils, 
respectively; 0.75- mm diam. wire). Bi: oversized induction 
coils (75xl25-mm size, 4-mm horizontal width, 52.5-mm av- 
erage spacing, 25 turns each; 0.4-mm diam. wire). PU: pickup 
coils, and CPU: extra detection coils for rf noise compensa- 
tion - see text - (10x18- mm diam., 6- mm width, 440 turns 
each; 0.2- mm diam. wire). A thin filling tube is attached to 
the aerogel container (here depicted with cylindrical shape). 

rounded by a closed copper box (23x16x10 cm 3 , thick- 
ness of 2 mm, not shown) for passive radio- frequency (rf) 
shielding. Two untuned, oversized rectangular induction 
coils are used to generate a very uniform rf field over 
the sample. Typically, a 90° tip angle is achieved by a 
0.5-ms rf pulse with 11.5-V rms excitation. 180° pulses 
induce negligible magnetization losses, less than 10% af- 
ter 1000 echoes refocused without applied gradient. Two 
circular detection coils lie in conventional Helmholtz con- 
figuration close to the sample ("PU" coils in Fig. 1). We 
use a pair of identical coils ("CPU" coils in Fig. 1), also 
implemented next to the sample and connected in op- 
position with the previous ones, to actively shield the 
device from pick-up noise. These two extra coils de- 
tect a noise signal from distant rf sources that almost 
exactly compensates for that induced in the main detec- 
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tion coils. Induction and detection coils are oriented so as 
to minimize the rf cross-talk, i.e., mutual induction. The 
set of detection coils is tuned to Larmor frequency with 
a parallel capacitor, and the quality factor is moderate 
(Q=43). The measured sensitivity corresponds to an ex- 
pected 90° FID signal of 43-^V rms amplitude per mbar 
pressure and cm 3 of 3 He gas with 100% nuclear polariza- 
tion. Two pairs of rectangular coils, connected in series, 
are designed to produce a uniform horizontal magnetic 
field gradient (25 mT/m/A) parallel to the static field 
axis. Computed relative variations over 3 cm around the 
center do not exceed 3%. 

We use a dual-phase lock-in amplifier (LIA; model 
7220, PerkinElmer Instuments, Oak Ridge Tennessee, 
USA) both to collect the NMR signal and to generate 
the resonant rf voltage. Rf pulse gating, dcphasing, and 
amplification are handled by a personal computer (PC)- 
controlled home-built circuitry. Gradient and rf pulse 
sequencing, as well as data collection via an A/D con- 
verter (sampling frequency: 10 kHz), are managed using 
a laboratory-written software. 



C. 3 He hyperpolarization system 

3 He gas is hyperpolarized inhouse using the metasta- 
bility exchange optical pumping technique. We use a 
home-built system that involves the basic units indicated 
by the block diagram in Fig. 2. The gas is polarized at a 
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FIG. 2: Block diagram of the hyperpolarization system. L: 
laser source. GHS: gas handling system. PL: pump laser 
beam. OP cell: optical pumping glass cell. C: peristaltic 
compressor. N: Nitrogen gas supply. V: vacuum system. PG: 
pressure gauge. NMR SB: shielding box surrounding the sam- 
ple and NMR setup sketched in Fig. 1. 

pressure of about 2 mbars in a cylindrical glass cell (50- 
cm length, 6-cm diameter). The 1083-nm light source 
is a 2-W ytterbium fiber laser (model OI-FL1083-2W, 
Keopsys, Lannion, France). The OP beam is expanded 
and collimated (2-cm full width at half maximum), and 
circularly polarized. A mirror back-reflects it to enhance 
light absorption by the small fraction (typically 1 ppm) 
of 3 He atoms excited to the metastable state by a weak rf 
discharge (3-MHz frequency, 2-W power). A dedicated 
peristaltic compressor is used to extract the polarized 
gas from the OP cell and transfer it into the sample con- 
tainer [47]. Upstream of the OP cell, the gas handling 
unit (not shown) includes a 3 He bottle (purity grade 
higher than 99.99%), a getter to keep impurities level 
well below lppm (model GC50, SAES Getters, Lainate, 



Italy), a high vacuum system (a turbomolecular pump 
backed by a membrane pump, electropolished stainless 
steel tubes and valves, pressure transducers), and an all- 
metal mass-flow regulator to control the gas flow into 
the OP cell. Downstream of the OP cell, the polarized 
gas is in contact with non-magnetic materials only (glass, 
Teflon, low-degassing plastic tubes, etc.). Leak-tight con- 
nections with respect to room air are required to avoid 
02-induced relaxation [48]. Dead volumes are kept as 
small as possible, with the constraint that excess relax- 
ation can result from drastic enhancement of surface-to- 
volume ratios. 

A weak static magnetic field (2 mT) is generated by a 
set of six square coils of horizontal axis, designed to pro- 
vide an optimal field homogeneity both for OP and NMR 
in a compact setup (Fig. 3). Relative inhomogeneities of 
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FIG. 3: Sketch of the He hyperpolarization system (front and 
side views), showing the arrangement of OP cell, compressor 
(C), driving motor (M), and of the six square coils (aver- 
age side length 40 cm, section 1.6x1.6 cm 2 ) that provide the 
static field (0.8516 mT/A at the centre). Field map is opti- 
mized both for OP and NMR by coil positions (±5.9, ±19.1, 
and ±37.8 cm) and number of turns (85, 100, and 224, respec- 
tively). 2x30 additional turns are wired in opposite directions 
onto the outer square coils, yielding a 1.3 mT/cm/A shim- 
ming gradient over the sample. A mirror (m) back-reflects 
the laser beam (lb) into the OP cell, open at both ends for 
gas inlet (i) and outlet (o) . 



the applied field do not exceed a few parts per thousand 
over the entire cell volume for OP, and a few ppm over a 
7-cm diameter spherical volume for NMR. The stray field 
inhomogeneity induced by the laboratory environment is 
weak, and local compensation of the dominant residual 
gradient is achieved using two shimming coils. 

Absolute measurements of the nuclear polarization in 
the OP cell are performed by optical polarization analy- 
sis of the fluorescence light emitted by the plasma. With 
2-W laser power, nuclear polarizations up to 60% can 
be obtained with no gas flow, and within the 30%— 40% 
range around 1 seem (standard cubic centimeter per 
minutc)(0.745 /Ltmol/s) flow rates, mostly due to pres- 
sure increase in the OP cell. The polarization losses in 
the compressor and connecting tubes, characterized in in- 
dependent measurements [49], typically amount to 10%, 
at best. 
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D. Probe gas management 

The filling tube attached to the aerogel sample is con- 
nected to a Luer polycarbonate stopcock valve located 
just outside the NMR shielding box. There, pressure (rel- 
ative to the atmospheric pressure) is measured with a 1% 
accuracy in the 1 mbar— 1 bar range using a piezoresistive 
sensor mounted in a nonmagnetic (plastic) case (Model 
24PCCFA2G, Honeywell, Morriston NJ, USA). For low- 
pressure diffusion measurements, pure 3 He is used. The 
gas flows directly from the OP cell and is accumulated 
inside the aerogel by the compressor. For high-pressure 
diffusion measurements, accumulation of large quantities 
of polarized gas is of limited interest considering the bot- 
tleneck set to magnetization buildup by wall relaxation 
(see Section IIIB). Therefore, N 2 is used as a neutral, 
nonrelaxing buffer gas to slow down the diffusion of the 
3 He atoms in the void spaces of the aerogels. The 3 He 
partial pressures vary from 10 to 200 mbars in the mea- 
surements, but most frequently no more than 30 mbars 
of polarized gas are introduced in the sample contain- 
ers (i.e., only a few //moles). Total pressures range from 
10 mbars to 1 bar. 

Before the NMR acquisition, the sample is flushed with 
pure N 2 (99.995% grade), then evacuated. Filling with 
polarized 3 He using the peristaltic compressor takes a 
few tens of seconds. Pure N 2 is optionally added to reach 
the targeted pressure value, then the measurement is per- 
formed. Right after signal acquisition, the gas is removed 
and the aerogel is kept under vacuum until the next mea- 
surement. Each series of measurements can last from a 
few hours to several days. When it is completed, the 
sample is open to room air and stored at ambient pres- 
sure. 



III. NMR MEASUREMENTS 

A. FID signals and line shapes 

The performed measurements do not require an abso- 
lute calibration of the NMR detection system. However, 
we have used a small transmitting rf coil located at the 
center of the setup to probe the global pick-up sensitiv- 
ity, and also to check that the 3 He signals correspond to 
nuclear polarization rates in agreement with those pro- 
duced by the optical pumping system. The measured 
noise level is about 0.7-//V rms when the LI A time con- 
stant is set to minimum (160 /lis). Nuclear polarizations 
vary from a few percents to 20% at best, due to relaxation 
in the aerogel, and the best spin-echo signal-to-noise ra- 
tios (SNRs) achieved are on the order of 30—50. Higher 
polarizations might be obtained if polarized gas is accu- 
mulated in a small glass cell, then rapidly expanded into 
the sample container. However, the procedure used here 
seems more suited to systematic and repeated measure- 
ments in the fragile silica substrates, and it still yields 
satisfactory accuracy. 



In our experiments, FID signal lifetime is limited by 
residual field inhomogencitics due to the magnetic envi- 
ronment of the experiment after shimming (Section II C). 
The motional narrowing regime can then be fully reached 
for low pressure gas, and exponential decay is observed 
with a time constant T 2 * exceeding 2.5 s at 60-mbar pres- 
sure (whereas T 2 * <0.5 s without shimming). At 1-bar 
pressure (no motional averaging) signal lifetime remains 
shorter, typically on the order of 0.2 s. 

For the diffusion measurements, a strong gradient is 
applied so that at all pressures the FID lifetimes be- 
come much shorter and the line shapes, obtained by fast 
Fourier transform (FFT) of the NMR signals, are those 
imposed by the one-dimensional (ID) spatial distribu- 
tion of resonance frequencies (see Section IIIC). Oth- 
erwise, line distortions occur due to diffusive relaxation 
or spectral edge enhancement effects [50, 51] and a va- 
riety of spectral profiles can be observed, depending on 
gradient amplitude (i.e., frequency span over the sample 
size) and gas pressure (i.e., diffusion time during obser- 
vation) conditions. Particularly important in highly dif- 
fusive systems, these line distortion effects have recently 
been systematically characterized with good SNR using 
hyperpolarized noble gases [31, 52]. Care has thus been 
taken to select gradient strengths and echo times so as 
to avoid systematic biases in the diffusion measurements 
[50, 53-56]. 



B. Longitudinal relaxation rates 

The contribution of dipole-dipole relaxation to the lon- 
gitudinal relaxation time T\ is totally negligible for pure 
3 He (relaxation time inversely proportional to the 3 He 
pressure, exceeding 800 h at 1 bar [57]). Only severe 
contamination by air may play a significant role in our ex- 
periments, through strong dipole-dipole relaxation by the 
paramagnetic 2 molecules (12.5 s at 1 bar and 21% oxy- 
gen concentration [48]). Bulk relaxation due to atomic 
diffusion in our residual field gradients is also negligible. 
Wall relaxation is thus the dominant process, and the 
magnetization lifetime depends on the nature of the ma- 
terials in contact with the gas, and on the area of exposed 
surfaces. 

The longitudinal relaxation time T\ is measured to be 
370 s in the empty machined PMMA cell. This corre- 
sponds to a specific relaxation time 6 S = T\A/V (A: 
area, V: volume) equal to 0.19 x 10 6 s/m, slightly lower 
than typical values reported for Pyrex glass (6 S — 0.5 — 
3.8 xlO 6 s/m [58-60]) and quartz (0 S = 0.3-2.8xl0 6 s/m 
[58, 61]). 

Surprisingly long T\ are measured in the aerogels, es- 
pecially in the commercial ones, despite their huge nom- 
inal specific areas (a few tens of m 2 per cm 3 ) and un- 
controlled cleanliness. T\ reaches 130 s and 200 s in 
samples A and J 3 , respectively. The contribution of 
the container walls can be estimated from the refer- 
ence PMMA cell T\. This yields particularly weak spe- 
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cine relaxation times on the silica aggregates (9 S = 6.6 
and 25 x 10 9 s/m, respectively), two orders of magni- 
tude longer than the record values obtained with special 
aluminosilicate glasses (low-boron content Corning 1720 
glass: 2.0xl0 7 s/m [62]; Schott iron-free Supremax glass: 
1.6-4.2xl0 7 s/m [63, 64]). T x does not exceed 26 s in 
the custom-made sample M, which is less than initially 
measured at 2.3 mT (58 s) but more than subsequently 
obtained at 0.1 T (19 s) in the same sample [37]. Fer- 
romagnetic impurities in the material [26, 58, 60] may 
be responsible for these hysteretic 7\ changes, also ob- 
served recently in bare glass cells exposed to high mag- 
netic fields [64] (a phenomenon reported for rubidium- 
coated glass cells only, so far [65, 66]). Using an esti- 
mate for the Pyrex T\ , and an approximate specific area 
(22.9 m 2 /cm 3 for a 98.2% similar custom-made aerogel 
[17]), the 58-s T\ initially measured in sample M would 
still yield 9 S = 1.5 x 10 9 s/m for the silica aggregates 
(shorter than, but comparable to, that of the commer- 
cial aerogels). In all samples, magnetization decays are 
indeed slow enough to allow the one-shot CPMG diffu- 
sion measurements described hereafter. 



C. Spin-echo attenuation 

A 90° rf pulse (0.5-ms duration) is followed by a train 
of 180° rf pulses (1-ms duration) to refocus the transverse 
magnetization at a regular interval Top- CPMG mea- 
surements are made for different values of Tqp (6—31 ms) 
and of applied gradient G (0—2 mT/m). Spatially non- 
selective tipping is obtained with a good angle accuracy 
by switching off (on) the applied gradient 300 [is before 
(after) the beginning (the end) of the rf pulses. 

For the CPMG sequence, the transverse magnetization 
has a decay time that can be written as 



2 1 2 J2,diff 

The first term on the right-hand side is the intrinsic re- 
laxation rate of the system, that for the gas simply cor- 
responds to the longitudinal decay rate (T2 — T\). This 
has been checked in the absence of applied gradient, using 
spin-echo trains with very short Top. The second term 
represents additional relaxation resulting from diffusion 
in a nonuniform magnetic field, due to an applied field 
gradient and/or the residual magnetic inhomogeneities 
after shimming. 

For a uniform applied gradient G, T^diff can be ex- 
pressed in the case of free gas diffusion as 

T2 ' diff = D^kGTcvf (2) 

where D is the coefficient of diffusion and the factor k 
corrects for the time variation of the applied gradient, 
which is switched off during the rf pulses. For a rectan- 



gular gradient pulse, this correction factor is: 

k = yV(3T CP - 26)/T* p (3) 

where 5 < T^p is the duration of the gradient [67] . 

All diffusion-weighted spin-echo trains are observed to 
decay exponentially An example of spin-echo train and 
FFT is presented in Fig. 4. values are obtained from 
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FIG. 4: a: Example of spin-echo train obtained in aero- 
gel A (with 33.4 mbars of 3 He plus 616.5 mbars of N2, 
G = 0.314 mT/m, and T C p = 21 ms). Saturation of the de- 
tection system by the strong rf pulses occurs between echoes 
(data points removed). Inset: Blow-up plot of first and second 
echoes, b: Fast Fourier transforms (FFTs) of the consecutive 
echo data sets, identified by individual symbol shapes, c: Ex- 
ponential fit of FFT decay, yielding T' 2 = 0.29 s. 

monoexponential fits to the square of the echo ampli- 
tudes, or the square of their Fourier components. The use 
of power data, rather than magnitudes, yields a more re- 
liable result as they are less sensitive to biases introduced 
by noise base lines [68]. 

The diffusion measurements are typically performed 
with short repetition times (Tcp= 6 or 11 ms) and fi- 
nite gradients (G=0.4— 1.2 mT/m), and ranges from 
25 to 1500 ms in aerogel (25 to 86 ms in the refer- 
ence cell). The contribution of intrinsic T2 to trans- 
verse relaxation is thus negligible at all pressures, and 
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T' 2 = T2.diff- The frequency spectra of the recorded 
echoes provide ID images of the gas magnetization in 
the samples (Fig. 5a). The shape of sample A is fairly 
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FIG. 5: Pixel-per-pixel analysis of frequency spectra, through 
monoexponential fits of the FFT decays, for spin-echo trains 
obtained in aerogels A (open symbols) and J2 (solid symbols) 
with G=1.56 and 0.43 mT/m, respectively, and Tcp=11 ms. 
a: Local initial amplitudes (along gradient axis). Statistical 
error bars, not plotted, are less than 2% of the maximum 
amplitude, b: Local decay times (left axis: A; right axis: J2). 
Statistical error bars are smaller than symbol sizes, except 
where plotted (i.e., on sample edges). 

well reproduced considering the selected spatial NMR 
resolution 27r/(7GT obs )=2.5 mm ( 7 =27rx32.43 MHz/T 
is the gyromagnctic ratio of 3 He nuclei, and T b s the 
observation time that is very close to the gradient du- 
ration 6). Resolution is lower for the displayed J2 profile 
(7.8 mm), but typical of most measurements presented 
in Section IV. The frequency spectra of the echoes also 
provide ID maps of the decay time (Fig. 5b). Spatial 
variations (±10%) are small enough to accurately yield 
a single exponential T 2 f° r the average magnetization. 

From each measurement, a diffusion coefficient D is 
extracted using Eq. (2). Experimentally, the I/T2 data 
can be checked to exhibit the expected linear variation 
with (kGTcp) 2 at all pressures above 100 mbars. Fig. 6 
displays, for instance, two series of measurements per- 
formed in the reference cell and in aerogel M at fixed gas 
pressures. Various combinations of parameters G and 
Tqp are used here, and in both cases all the data col- 
lapse on a straight line. In the empty reference cell, 
D coincides with the free-diffusion coefficient -D gas of 
the 3 He atoms in the probe gas. In aerogel, D corre- 
sponds to an effective diffusion coefficient D e ff that de- 
scribes the restricted diffusion of the atoms within the 
porous medium. Below 100 mbars, moderate departures 
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FIG. 6: Compilation of spin-echo decay rates measured in 
the reference cell (solid squares) and in aerogel M (open cir- 
cles) at fixed gas pressures. Diffusion-weighting parameter 
(kGTcp) is varied through changes of gradient strengths (G: 
0.07-2.2 mT/m) and echo times (T C p: 4.4-31 ms). The ex- 
pected linear dependence [Eq. (2)], observed over more than 
two decades, yields D gas = 1.148 ± 0.004 cm 2 /s in the refer- 
ence cell, and D c g — 0.595 ± 0.008 cm 2 /s in aerogel M, for 
the two 3 He-N2 mixtures used (P^ = 1.41 bar and 1.70 bar, 
respectively - see Section IV A). 



from the (kGTcp) 2 dependence can be observed in aero- 
gel at high GTcp, corresponding to reduced decay rates 
l/Tj (hence, overestimated values of D), typically above 
kGT C p = 4 x 10~ 6 sT/m as reported in Rcf. [39]. 



IV. RESULTS 



A. Free gas diffusion 



The free-diffusion coefficient D gas for a 3 He gas with 
partial pressure P3 in a mixture containing N2 with par- 
tial pressure Pn2 can be written as: 



N2 



V Z {T) V m2 {TY 



(4) 



where the reduced diffusion coefficients can be expressed 



as X> 3 = 1.967 x 

vl.65 



(T/300) 1 71 and V m2 



0.811 x 

(T/300) 1 ' 00 in units of atm x cm 2 /s (1 atm = 1.013 x 
10 5 Pa), with the temperature T in Kelvin [36, 69]. 
Therefore, the scattering on N 2 molecules leads to an ef- 
fect of the gas composition on atomic diffusion that can 
be accounted for using the equivalent pressure of pure 
3 He gas: 



P3 + P 



N2 



V 3 (T) 
V m2 (T) : 



(5) 



with P3/D3N2 = 3.462 at 293 K. This allows to merge 
the results obtained in 3 He— N2 mixtures and in pure 3 He 
gas. For free diffusion, 1/-D gas and hence T' 2 are expected 
to linearly depend on P^ q . 

CPMG measurements performed in the reference cell 
are presented in Fig. 7. These data are obtained at 
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various 3 He pressures, with several combinations of pa- 
rameters G and Tqp- They actually all collapse on a 
straight line when plotted as a function of the equiva- 
lent 3 He pressure P^ q . Moreover, the measured slope 
5.06 ± 0.04 mbar -1 x s/m 2 (statistical error bar) is 
in fair agreement with the expected value (l/£>3 = 
5.23 mbar -1 x s/m 2 ) given the 4% uncertainty due to 
long-term fluctuations of atmospheric pressure and tem- 
perature. This provides a consistency check for our appa- 
ratus, gradient assignment, and measurement technique. 



1- 



E 

^ 0.1 
XL 10 



4 . 

O >< 



3* 



100 



1000 



10. 




1000 2000 3000 
P„ eq (mbar) 



FIG. 7: Comparison of diffusion measurements performed in 
free gas, and in aerogels A and M. The same three sets of 
reduced data 1/D e ff, plotted against pressure P^ 1 , are dis- 
played in linear and logarithmic scales. Stars: reference cell 
data, free diffusion. Circles (respectively, squares): aerogel M 
(resp. A) data. Crossed circles: previous measurements in 
aerogel M [38, 39]. Solid line: linear fit of the free-diffusion 
data (slope: 5.06±0.04 mbar -1 x s/m 2 ), in good agreement 
with the expected P^/'D^ dependence. Dotted lines: phe- 
nomenological power-law fits with exponents 0.47 (aerogel A) 
and 0.60 (aerogel M). 



B. Gas diffusion in aerogel 

Results obtained in aerogel samples A and M are also 
displayed in Fig. 7, so that they can be directly compared 
with the free-diffusion results. For clarity, the results 
obtained in samples Ji— J 4 are presented separately in 
Fig. 8. 
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FIG. 8: Reduced diffusion data 1/D c g as a function of Pg q 
for aerogel samples Ji (plus signs), J 2 (open triangles), J 3 
(open diamonds), and J4 (crosses). Comparison is made with 
the linear fit of free-diffusion data (straight line), and power- 
law fits of aerogels A and M data from Fig. 7 (dotted lines). 
The Ji— J4 data can be consistently fitted by a power-law 
functions with exponents 0.53 or 0.64 (see Table I). 



In Fig. 7, data from earlier measurements in sample M 
are also included for comparison. These measurements 
[38, 39] have been performed in a 0.1 T whole-body im- 
ager, using the broadband rf channel and the gradient 
system of that clinical unit. Homemade antennas were 
used for induction and detection at 3.24 MHz, as well 
as the in-house 3 He-NMR software for sequence man- 
agement and data acquisition developed for in vivo lung 
studies [35, 70]. Pure 3 He gas was used for these CPMG 
measurements (G = 0.21 mT/m, Tqp = 10 -2 s). De- 
spite all these differences in apparatus and experimental 
conditions, the two sets of data are perfectly consistent. 
The present work significantly extends the probed range 
of gas pressures. It confirms both the global decrease in 
diffusion coefficient (as compared to free diffusion) and 
the strong nonlinear dependence on gas pressure. These 
results therefore appear to be quite robust, and reliably 
associated to intrinsic bulk characteristics of the aerogel 
(on which time or rough pressure variations seem to have 
no effect). 

In Fig. 7, the results obtained in the other aerogel (A) 
are qualitatively similar but show a strikingly stronger 
departure from the free-diffusion behavior, revealing a 
much more restricted 3 He diffusion. This may be at- 
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tributed to the difference in aerogel porosities, the silica 
volume density being twice larger in A than in M. The 
quantitative difference in pressure dependences is empha- 
sized by the logarithmic plot, that shows that the two sets 
of data are best fitted with distinct power laws. 

Type J aerogel samples have exactly the same nominal 
porosity as A, and yet the data for Ji— J4 are actually 
much closer to those of M (see Fig. 8). Finally, data ob- 
tained in the four samples of type J are quite satisfacto- 
rily comparable, revealing only minor sample-to-sample 
variations (again, there is no sign of individual alteration 
by inflicted mechanical stress). These data are of lower 
quality, the SNR being lesser due to the smaller sizes of 
the Ji— J4 samples. 

1 / Z?eff exhibits a strong nonlinear pressure dependence 
in all aerogels. The logarithmic graphs in Figs. 7 and 8 
show that the data can be described by simple power- 
law functions, with exponents lying in the 0.48—0.65 
range. For aerogels A and Ji— J4, the data depart from 
that common behavior (longer are measured) below 
100 mbars, and this value is chosen as a lower bound for 
the fitting range. At the lowest pressures, limits of the 
validity domain of Eq. (2) (see Section III C) are indeed 
reached, due to the experimental conditions for measure- 
ments performed in these samples (10 6 /cGTcp = 4.5 and 
6.1 Ts/m, respectively). With the choice made for sam- 
ple M (10 6 fcGXcp = 3.0 Ts/m), no such departure occurs 
down to 10 mbars and all data are included in the power- 
law fit. Fit parameters are listed in Table I to sum up 
samplc-to-sample variations. 



V. DISCUSSION 

The use of hyperpolarized 3 He is demonstrated to al- 
low diffusion measurements with good accuracy, over a 
wide range of gas densities. The compact polarizer is 
perfectly suited to online preparation of the small gas 
quantities required to test the aerogel samples, and sys- 
tematic measurements are easy to carry out and re- 
peat. The high porosity samples under study are se- 
lected among silica aerogels that are most frequently 
used as host porous matrices for phase transition stud- 
ies in low-temperature physics. None of these particular 
samples has been submitted to cool down- warm up, or 
compression-decompression, cycles comparable to those 
routinely performed in the low-temperature experiments. 
Exposure to high pressure or thermal stress has recently 
been demonstrated to induce irreversible changes in aero- 
gels, that may macroscopically affect their structure [22]. 
It might thus be interesting to perform nondestructive 
NMR diffusion measurements in aerogel both before and 
after use at low temperature. 

The consistent set of robust NMR data collected in the 
differently produced virgin samples reveals common bulk 
aerogel features, that are analyzed in terms of specific 
properties of the microscopic structure in this section. 
The unusually high porosities of the samples allow a par- 



ticularly simple description of the gas kinetics inside the 
aerogels. Restricted diffusion is more often studied in 
denser media, where atoms or molecules are confined in- 
side an ensemble of closed or interconnected pores that 
strongly constrains, by its geometry or tortuosity, the 
thermal motion of the probe species. In this work, the 
solid matrix occupies a very small fraction (1.5% to 3 %) 
of the volume of the aerogel, and the hyperpolarized 3 He 
atoms diffuse inside the widely accessible, huge space left 
open by the sparsely distributed silica aggregates. 



Decay times of spin-echo amplitudes are measured to 
be two to five times longer in aerogel than in free gas. 
The decrease of the atomic diffusion coefficient is quite 
significant at pressures exceeding 1 bar, indicating that 
atom-atom and atom-silica collisions have comparable 
frequencies at these high gas densities. Surprisingly, a 
strong pressure dependence of D e g is still observed at 
low pressures. As the probability of atom-atom scatter- 
ing vanishes in a very dilute gas, the relative effect of 
pressure variations on the 3 He diffusion coefficient would 
be expected to become increasingly small at low pressures 
(and some characteristic length scale, if any, of the silica 
network should determine the finite value approached by 
the atom mfp at zero pressure - see Section VA). It is 
yet measured to remain large at low pressures, indicating 
that the free motion of the atoms is indeed allowed over 
large distances in aerogel (at 10 mbars, for instance, the 
mfp for free diffusion in 3 He is on the order of 40 /xm 
[71]). This is systematically observed, in the custom- 
made aerogel sample and in the set of commercial ones 
that are obtained from two distinct manufacturers. It 
brings further support to the suggested nonuniform dis- 
tribution of density of silica strands inside the aerogels, 
mentioned in the Introduction, that had been put for- 
ward to explain some preliminary results [38, 39]. 



The NMR gas diffusion data actually have two distinct 
features that compel description of the aerogel by more 
than a single characteristic length: the dominant power- 
law increase with pressure, and the departure from it 
that can be noticed at very low pressures. This depar- 
ture has already been discussed in Ref. [39], and shown 
to be enhanced at large fcGTcp. As pressure decreases, 
the atomic diffusion gets more and more sensitive to scat- 
tering by the silica strands. This enhances the effects of 
spatial distribution inhomogeneities and the multiexpo- 
nential nature of the spin-echo decay, as in other porous 
media [25]. The longest lived mode is favored for all 
length scales below the gradient-induced resolution, and 
monoexponential fits systematically lead to notably un- 
derestimated relaxation rates (i.e., increased 1/Z? e ff val- 
ues). The following discussion focuses on the power-law 
dependence on pressure observed in all samples and on 
its analysis. Results are also compared to those obtained 
with other investigation techniques. 
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Sample 


Manufacturer 


Porosity 
[%] 


Specific area 

r 2/1 

K/g] 


Exponent 
j3 (stat. err.) 


Prefactor [value for P.^ —1 bar] 
a (stat. err.) [s/cm 2 ] 


M 


custom made 


98.5 


- 


0.60 (0.01) 


1.25 (0.01) 


A 


Airglass 


97 


465 


0.47 (0.01) 


3.16 (0.03) 


Ji 


Matsushita 


97 


510 


0.53 (0.02) 


1.42 (0.02) 


J 2 


Matsushita 


97 


510 


0.53 (0.02) 


1.54 (0.02) 


Js 


Matsushita 


97 


510 


0.64 (0.01) 


1.31 (0.01) 


J 4 


Matsushita 


97 


510 


0.64 (0.02) 


1.27 (0.02) 



TABLE I: Characteristics and fit parameters for the six aerogel samples. Exponents (/3) and prefactors (a) are obtained by 
phenomenological fits of the NMR diffusion data to 1/D c s = a x (P^) 13 , where \/D e g is in s/cm 2 and P^ q in bars (see text). 



A. Crude model: uniform medium with a single 
characteristic mfp 

The 3 He diffusion process leading to spin-echo attenua- 
tion in the applied magnetic-field gradient can be viewed 
as a random walk in which free-flying atoms take suc- 
cessive individual steps that are terminated either by a 
collision with another gas atom or molecule ( 3 He or N 2 ) 
or with silica. The frequencies of the two types of col- 
lisions being additive, the mfp A c ff of this random walk 
can be related to the mfp A gas for gas collisions and A aero 
for collisions with the silica strands in a straightforward 
manner: 



1 

A e ff 



1 



A. 



+ 



gas 



(6) 



The spatial distribution of silica is here assumed to be 
homogeneous, and the global effect of the aerogel struc- 
ture on atomic motion is thus described by the single 
parameter A aoro . The effective diffusion coefficient D c s 
can then be written as: 



1 



1 



+ 



1 



(7) 



where the first contribution (gas-gas scattering) involves 
the free-diffusion coefficient of 3 He introduced in Sec- 
tion IV A, and the second one ( 3 He— silica scattering) 
introduces the constant diffusion coefficient Z? aoro = 
^A aoro /3, where v is the average thermal velocity (v = 
^SksT/ivm; ks- Boltzmann constant, m: atomic mass, 
T: temperature). The quantity l/D c g should therefore 
increase linearly with the gas pressure (or the equivalent 
pressure Pg q for a 3 He — N 2 mixture), starting from a 
finite zero-pressure value that may vary from sample to 
sample with the degree of confinement of the 3 He atoms 
by the silica network. 

Over the wide range of the gas pressures investigated, 
the expected linear dependence on pressure is clearly not 
observed. This strongly rules out the description by a 
uniform density of silica scattering centres in the aero- 
gels and an average mfp. Instead, a broad distribution 
of mean free paths at scales smaller than the imaging 
resolution has to be assumed to account for the experi- 
mental observations. This is quantitatively discussed in 
the following (Section VB), where the main features of 



required mfp variations are correlated, over the probed 
range (0.1—40 /xm), with the measured pressure depen- 
dences. 



B. Phenomenological model: continuous 
distribution of mfp's in aerogel 

To describe diffusion in heterogeneously distributed sil- 
ica, a continuous distribution / of mfp's is introduced. 
The probability that the 3 He-silica mfp A aoro lies be- 
tween A a cro and A aoro + rfA acro in some volume fraction of 
the sample is /(A a ero)dA aoro . Assuming that the atoms 
rapidly exchange between regions of space characterized 
by different mfp's, the 3 He magnetization is uniform and 
decays with a weighted average of local relaxation rates: 



'2,diff 



/ (A a0 ro )^A ao 
?2,diff(A aero 



Since the local relaxation times are given by: 



£2,diff 



(-^aero) 
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(~fkGT( 



CP) 



1 



D„ 



+ 



1 



-^aero (^aero ) 



(8) 



), (9) 



the effective diffusion coefficient derived from the global 
decay time T 2j diff using Eq. (2) reads: 



D, 



off 



/ 



f (^aero ) -^aero ( ^aero ) -^gas ^-^acro 



H - -^aero 

/(A 



(•^aero) 

)T> 3 d\ a 



-faero(^aero) H" -£3 



(10) 



The physical quantity -P a oro(A aC ro) = 32?3/vA ae ro7 that 
appears in the right hand side of Eq. (10), has the di- 
mension of a pressure. It represents the contribution of 
the 3 He-silica scattering process to hindering of atomic 
motion, on the A aor o length scale, by the solid matrix. 
Eq. (10) shows that the integrand is a simple rational 
function of the probe gas pressure. The variations of 
D c s with Pg q therefore put strong constraints on the 
mfp distribution / introduced in this model. 

Experimentally, D e s is proportional to (P^)^ 13 over a 
wide range of pressures. This is mathematically equiv- 
alent to set D e g(aP^ q ) — a - * 3 D c ff(P^ q ), a being is an 
arbitrary scaling constant. Through a straightforward 
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change from A aer o to aA acro as dummy variable in the in- 
tegral [Eq. (10)], this relation can be shown to hold for all 
values of P 3 eq if, and only if, /(A acro /a) = a 2_/3 /(A aoro ). 
This translates into /(A aoro ) oc (Xacro) 13 ^ 2 - The mfp dis- 
tribution, itself, corresponds to a power law, and its expo- 
nent is directly related to that measuring the 1 / D e ff data 
scaling rate versus pressure. Although P^ q is varied over 
more than two decades, D c ff is actually measured over 
a finite pressure interval. Additional information would 
thus be required to fully characterize /(A aoro ). Yet, the 
distribution /(A aoro ) oc (A aoro ) ,3 ~ 2 provides a convenient 
tool for quantitative analysis and discussion of the exper- 
imental results. 

The effective diffusion coefficient can easily be com- 
puted from Eq. (10) using the distribution function 
/(A aoro ) (A a0 ro) _2+,3 1 provided that / is suitably nor- 
malized. Since exponents typically range from 0.5 to 
0.7 in the aerogel samples under study, no divergence 
occurs at infinitely large A, but some cutoff must be in- 
troduced at small A to keep J f(X)d\ finite. Numerical 
calculations are performed to probe the influence of both 
the upper and lower bounds of the mfp distribution, and 
to examine the constraints experimentally set by _D er r and 
its variation with P^ q in the 10 mbar— 3 bar range. For 
simplicity, /(A aoro ) is assumed to be zero outside the in- 
terval I < A a ero < L. Fig. 9 displays the results obtained 
at fixed (i, arbitrarily set equal to 0.5, for different val- 
ues of parameters I and L. As expected, the broader the 
mfp distribution, the closer to a power-law dependence is 
the variation of -D er r with pressure. Variations of the up- 
per bound L lead to small changes of D e g that can only 
be noticed at low pressures. These changes are found 
to be quite negligible above 10 mbars for L greater than 
1 mm. The lower bound I has a much stronger impact. 
1/DeS increases when I is reduced, which corresponds to 
the aerogel hindering gas diffusion and restricting free 
paths to very short distances. On the contrary if, for in- 
stance, I is raised to 1 /an the atoms move almost freely 
inside the aerogel over the probed pressure range, and 
D e g approaches Z? gas . As the exponent of the pressure 
dependence exhibits no noticeable drift towards 1 at high 
pressure (log plot in Fig. 9), the experimental data def- 
initely rule out I > 30 nm: mfp's on the order on a few 
nanometers must be encountered in all the aerogels under 
study. 

Taking advantage of the weak impact of the upper 
bound L (all the more negligible since P^ q <100-mbar 
data have been discarded in the first of all samples but 
aerogel A), one can actually extract from the prefactor 
values a in Table I estimates for the smallest mfp encoun- 
tered in each sample. Computing 1 / D e g as a function of 
I using Eq. (10) and the fit exponents (3 (L being conser- 
vatively set to 3000 ^m), at fixed pressure P| q = 1 bar, 
adjustments to the fit prefactors a yield results with quite 
reasonable orders of magnitude: I —29, 24, 22 and 24 nm 
for samples Ji, J2, J3, and J4, respectively; I =31 nm 
for aerogel M; I = 6.5 nm only for aerogel A, where 3 He 
atoms experience the most restricted diffusion. 
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FIG. 9: Pressure dependence of 1/D e g computed from 
Eq. (10) for broad power- law mfp distributions in aerogel 
/(A aero 

) a (A aero) with f3 — 0.5, as a function of lower 
(I) or upper (L) bounds. The mfp distributions range either 
from / =0.01 fim to L=l, 3, 10, and 30 mm (symbols), or 
from I = 0.01, 0.03, 0.1, 0.3, and 1 fim to L = 30 mm (lines), 
as indicated in the legend. Dashed line: free diffusion limit. 
Inset: Blow-up of the log plot, focused on the small varia- 
tions induced at low pressures by changes in upper bound L 
(at constant lower bound I =0.01 /im). 



More numerical results can be found in the supplemen- 
tal material provided with this article through EPAPS 
[72] . In particular, shunt effects associated with the pres- 
ence of macroscopic aerogel-free spaces inside the exper- 
imental cells are quantitatively investigated. Their im- 
pact is small for volume fractions below 10%, and would 
mainly result in a downwards curve bending at low P^ q 
that is not observed experimentally. Other types of broad 
mfp distributions (e.g., Gaussian, bimodal, and exponen- 
tial distributions) are also easily ruled out, as computer 
calculations lead through Eq. (10) to a pressure depen- 
dence that is qualitatively very different from the ob- 
served one. 

In summary, the spin diffusion measurements per- 
formed in aerogel strongly suggest a description of this 
particular medium by a broad range of mfp's with power- 
law statistical distribution. The global effect of silica net- 
work and its complex structure is best modeled with 
/(A a oro) cc (A aoro ) _1/ , with 1.3 < v < 1.6 depending on 
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the investigated aerogel sample. Data analysis indicates 
that these distributions must extend at least from less 
than 30 nm to more than hundreds of micrometers. 



C. Comparison with other studies 

To our best knowledge, there is no comparable quan- 
titative study of gas (self-)diffusion in aerogels [73]. Re- 
stricted gas diffusion has been reported in a recent at- 
tempt to probe aerogel pore-structure with imaging and 
spectroscopy NMR techniques using thermally polarized 
129 Xe gas [74, 75]. A 3-fold decrease of the Xe diffu- 
sion coefficient is observed for pressures ranging from 1 
to 30 bars, in several small specimens of unspecified ori- 
gin and nominal characteristics. Other dense probe fluids 
can also be used to probe silica aerogels, e.g., water [76] 
or methanol [77]. But liquid phase NMR involves much 
shorter relaxation times and the quantitative analysis of 
the attenuation data is much less straightforward (spin- 
spin or spin-lattice relaxation processes must be taken 
into account [77]). Most reported 'gas diffusion' mea- 
surements [73] are gas phase studies that actually consist 
in transport measurements (permeametry), with assess- 
ment of gas flow through virgin [78-81] or partially den- 
sified [80, 82, 83] aerogels. In this case, a mean pore size 
of the interconnected open network is usually inferred 
from measured flow rates of gas through the pores using 
a suitable geometrical model of the open system. 

Flow measurements actually provide information 
about the structure of the porous matrix in the Knudscn 
regime, i.e., at low enough pressure, when free molecu- 
lar diffusion dominates. The gas-flow conductance then 
becomes pressure independent, and it is solely deter- 
mined by the cumulative effect of collisions on the silica 
walls [83-86]. The standard capillary-bundle model for 
porous media appears to give an acceptable description of 
both static or dynamic flow experiments in aerogel [78- 
80, 82]. Using a tortuosity coefficient t to account for 
the pore morphology and connectivity of the pore net- 
work, the diffusion coefficient is written as D e g = Dku/t 
with reference to the Knudscn diffusion coefficient Z?Kn in 
straight parallel, non-overlapping cylindrical pores of in- 
finite length. In this model D^n = vcj)/3, where the tube 
diameter cf) is equal to the characteristic size of the real 
pore structure under consideration (e.g., the hydraulic di- 
ameter, or mean intercept length). The <f> values derived 
from these molecular flow measurements range from 10 to 
50 nm, in agreement with the hydraulic diameters com- 
puted from standard stereology: = 4e/pE (e: porosity, 
p: density, and S: specific area) or the mfp's obtained 
from SAS measurements [80, 82]. And the effective dif- 
fusion lengths are 1.4 to 3 times larger than for straight 
line motion, due to tortuosity [79, 80]. 

Extending the capillary-bundle model to the aero- 
gels used in the present work, in the Knudsen regime 
one would expect diffusion coefficients on the order of 
0.6 cm 2 /s for samples of both type A and J, since the 



hydraulic diameters computed from their commercially 
specified specific areas and densities are quite compa- 
rable for these two (0=127 and 116 nm, respectively). 
The large differences in D c g experimentally observed by 
NMR at all pressures, in the two types of aerogel, sug- 
gest that either (1) the effective pore size derived from 
the capillary-bundle model suited for the analysis of flow 
measurements is not relevant to characterize the influ- 
ence of the complex ramified silica network on pure self- 
diffusion in the gas, or (2) the structure of the aerogel 
matrix in the high porosity samples under study is very 
different from that of most aerogels used in the above- 
quoted references. 

In several experiments, the Knudsen regime is reported 
to be reached for average gas pressures typically below 
1 bar for aerogel densities ranging from 0.1 to 2.2 g/cm 3 
[80, 82]. However, a clear transition from molecular to 
viscous flow has once been observed around 150 mbars 
in a 0.28-g/cm 3 sample (yet, the authors computed an 
effective pore radius on the order of 19 /jm, using a 
modified capillary-bundle model) [78]. Moreover, a sim- 
ilar observation can be found in one recent systematic 
study devoted to samples from an aerogel monolith grown 
with comparably high porosity (density 0.06 g/cm 3 ; area 
800 ±30 m 2 /g) [81]. Interestingly enough, the conclusion 
drawn from flow measurements performed in that study 
is that "the absence of any clear break (in the transition 
from diffusive to viscous flow) supports the assumption of 
a continuous distribution of pores sizes within the range 
of mean free paths in this experiment" . The experiment 
was actually performed with He mfp's ranging from 90 to 
3.5 /zm (measurements were also repeated with N 2 gas, 
for mfp's ranging between 30 and 1 /im, and a similar 
behavior has been observed). Furthermore, this anal- 
ysis correlated with an intriguing "significant pore vol- 
ume (detected) throughout the measuring range of the 
instrument from 2 to 200 nm" (by nitrogen desorption 
isotherm technique) [81]. Hence, the upper and lower 
boundaries of the continuous pore size distribution as- 
sumed for the high porosity aerogel specimens investi- 
gated in Ref. [81] are quite comparable to those derived 
in the present NMR work for 3 He (self-)diffusion mfp's, 
using the power-law distribution model presented in Sec- 
tion VB. 



VI. CONCLUSION AND PROSPECTS 

In summary, laser optical pumping has been combined 
with low field NMR to nondestructively probe high poros- 
ity aerogel samples, by spin-echo diffusion measurements 
with hyperpolarized 3 He gas. Analysis of the diffusion 
data has been based on a statistical description of the gas 
motion inside the solid matrix. Restricted diffusion has 
been related to 3 He mfp changes induced by interactions 
between the probe atoms and the sparsely distributed sil- 
ica strands, a gas kinetics approach particularly suited to 
the widely open aerogels under study. The experimental 
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findings are satisfactorily accounted for, assuming broad 
mfp's distributions inside the aerogels with ranges and 
power-law probability weights that are determined from 
the pressure dependence of the measured effective diffu- 
sion coefficients. 

The conclusions drawn from our study support the idea 
that these high porosity aerogels are not well-defined sys- 
tems. It may hence contribute to cast doubt on simple 
models and interpretations of low-temperature physics 
experiments, and fosters the growing concern that they 
may not be ideal host media when high control of in- 
duced confinement and disorder is required. Moreover, 
some aspects of the methodology and experimental tech- 
niques that have been used may be directly applied, or 
extended, to a variety of other porous systems, as follows. 

In our work, a buffer gas (N 2 ) has been used to slow 
down 3 Hc diffusion in the aerogel void spaces, and to 
achieve a significant increase of the pressure measure- 
ment range. This is applicable in all systems where the 
gas pressure cannot be increased ad libitum. For in- 
stance, diffusion imaging of inhaled hyperpolarized 3 He 
gas is currently used to observe in human lungs spatial 
variations of alveolar size related to physiology [36, 87], 
posture [88], respiratory conditions [89], or disease [89- 
91]. Operation at atmospheric pressure being mandatory 
in vivo, mixing with biologically inert heavy buffer gases 
may be advantageously used to enhance contrast for such 
studies (e.g., SF 6 decreases the 3 He diffusion coefficient 
by a factor up to 4 at 1 bar, depending on dilution) but 
also, more generally, to reduce diffusion-weighted attenu- 
ation by encoding gradients for improved SNR sensitivity 
and image quality [92]. 

Using hyperpolarized 3 He, we have been able to probe 
our aerogel samples over a very broad range of diffu- 
sion lengths with the same measurement technique. This 
provides a powerful alternative to standard NMR ap- 
proaches, often limited by technical aspects (e.g., acqui- 
sition rates, gradient amplitudes or ramping times), pro- 
vided that the pressure can be directly varied by several 
orders of magnitude in the porous medium. Applica- 
tion to denser silica aerogels is straightforward, but may 
be considered as well for other rigid porous systems al- 
lowing fast gas diffusion. For instance, diffusion mea- 
surements could be considerably extended in vitro, using 
fixed lungs with preserved air spaces [92], to bridge the 
gap between conflicting in vivo lung investigations per- 
formed on time scales ranging from a few milliseconds to 
tens of seconds (using gradient echo [89, 90], spin echo 
[35] and spin tagging [93, 94] techniques), and to clear up 
the link between 3 He NMR data and airways structural 
characteristics [95-97]. 

Laser polarization of the probe gas actually yields high 
SNRs to perform highly accurate NMR measurements in 
a variety of model systems (single pores [56, 98], rocks 
[99] or packed glass beads [99, 100]). The technique yet 
requires that a sufficiently large nuclear polarization is 
preserved inside the porous medium. This potentially 
limits the type of materials that can be probed, at least 



in the static mode operated here. However, for materials 
with moderately low specific relaxation times, continuous 
flow techniques can be used to maintain a suitable sta- 
tionary magnetization, through replenishment by freshly 
polarized gas. Most demonstrations involve hyperpolar- 
ized xenon gas so far [100, 101], but flow measurements 
are also very easy to perform with the compact 3 He gas 
production system described in Section II C. 

Finally, a specific interest of operation at low mag- 
netic field may be emphasized. Interpretation of NMR 
diffusion data is known to crucially depend on adequate 
background gradient compensation, but also on the pre- 
dominance of applied gradients over internal ones. In 
our work, local-field inhomogeneities due to the aerogel 
matrix were irrelevant for the analysis of the measure- 
ments. But NMR studies with hyperpolarized gas are 
highly sensitive, and they may equally reveal weak bulk 
magnetic susceptibility effects (e.g., the influence of glass 
containers at high field [31]) or artifacts generated by 
dominant susceptibility- induced field gradients (e.g., in 
ventilation imaging at 1.5 T [33, 34]). Operation at low 
or ultralow field can be very helpful to separate the effects 
of applied and local gradients for improved quantitative 
understanding in complex porous systems, as recently 
demonstrated in lungs [35, 36, 87]. 
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Supplemental material: EPAS Document 
E-JCPSA6-123-034528 [72] 

All the notations used in this supplemental material 
are those introduced in the related article. 

The first section investigates the effects of macroscopic 
void (aerogel-free) spaces present in the experimental vol- 
umes on the NMR diffusion measurements with hyperpo- 
larized 3 He gas. The probe atoms are strongly confined 
inside the aerogel, but diffuse freely (hence, much faster) 
in these void spaces. If the contribution of such diffusion 
"shunts" to the spin-echo attenuation is significant, the 
global effective coefficient D e ff extracted from the NMR 
data must differ from that of the bulk aerogel. Shunt ef- 
fects are quantitatively assessed, and discussed as a func- 
tion of the volume fraction of open void spaces. 

The second section provides examples of computed ef- 
fective diffusion coefficients in bulk aerogel, for a few 
kinds of broad mfp distributions. Results are successively 
presented for Gaussian, bimodal, and exponential distri- 
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butions. They contrast strongly with those obtained with 
the power-law distributions used in the article, that sat- 
isfactorily describe the experimental diffusion data. 

A. Shunt effects 

Experimentally, the void spaces open to 3 He diffusion 
may be large cracks that have developed inside the aero- 
gel samples, or gaps between the bulk material and the 
container walls for aerogel samples that do not fill up the 
whole cell volume. 

A finite fraction x of the sample volume where mfp 
Aacro gets infinite is introduced to describe such a macro- 
scopic "aerogel-free" region of space, where the atoms 
are free to move. As in Section V-B of the related arti- 
cle, fast exchange is assumed between this region and the 
bulk aerogel, leading to an effective diffusion coefficient 
that can be written as: 

D —(1 x) ^ •^(^ acro )^ 3 ^ acro (^ aor °) c ^ acro -^x^ 3 

J D 3 + -Daero( Aacro )^3 q P-^ 

(S-l) 

Fig. 10 shows the variations of 1/D e g with the void frac- 
tion x, for a fixed mfp distribution in aerogel (ranging 
from 10 to 30 mm, with (3 =0.5). For all pressures, triv- 
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FIG. 10: Computed pressure dependence of 1/D c ff as a func- 
tion of the void fraction x [Eq. (S-l)], for a mfp distribution 
in aerogel /(A aoro ) oc (Aacro)' 3-2 with j3 = 0.5, ranging from 
/ =10 nm to L =30 mm. 

ially, -D c ff increases towards -D gas as x increases. How- 
ever, large global changes of 1/D e s (i.e., of spin-echo de- 



cay times) occur only for x > 10%. The changes induced 
by the void fraction are relatively larger at lower pres- 
sures. At small x, the log plot indicates that they trans- 
late into a significant downwards bending of the curves at 
low Pg q that is not observed experimentally. At interme- 
diate x, the pressure dependence is strongly altered over 
the whole pressure range, and cannot be satisfactorily 
described by a single power-law variation. 

In the custom-made aerogel M (grown-in sample, at- 
tached to the glass walls) , the void fraction corresponding 
to the visually observed cracks is very likely to be ex- 
tremely small. Since no significant departure from the 
power-law fit is observed in the data down to P^ q — 
10 mbars, their actual contribution to D c ff is negligible. 

Void spaces are larger in the commercial aerogel sam- 
ples. Some have shape imperfections (e.g., small pieces of 
material missing on the corners, and irregular edge sur- 
faces) that leave macroscopic open spaces with dimen- 
sions below 1 mm. All samples, irreversibly shrunken 
after a few pressure cycles, no longer occupy the whole 
volume of the containers. Conservative estimates for the 
volume fraction x hardly reach 10%— 20%. The largest 
void fraction must be that of sample Ji (crushed sample 
with trapezoidal shape, already not matching its con- 
tainer in thickness). But, the exponent associated to the 
pressure dependence measured in aerogel Ji is actually 
smaller than those obtained in aerogels J2 and J4, and 
not larger. 

All this supports our belief that the measured features 
can trustfully be assigned to be those of bulk aerogel. 



B. Diffusion coefficients computed for various mfp 
distributions. 

The expected pressure dependence of 1 /-D e ff can easily 
be computed for any type of normalized mfp distribu- 
tion function using Eq. (10) of the related article. The 
expression of the effective diffusion coefficient may be 
slightly modified by simple algebra (using -P a ero(A a cro) = 
32? 3 /v Aacro), to let the rational function of pressure and 
mfp involved in the integrand appear more explicitly: 

= f /(A acro )P 3 rfA aero _ X> 3 f /(A) A dX 

(Aacro) + P 3 ° q Pi J Al + A (PP/Pl) 

where Pi is an arbitrary reference value for pressure, 
and Ai = 32? 3 /uPi the corresponding mfp for free dif- 
fusion (e.g., Ai = 399.2 ^m at 293 K for Pi = 1 mbar). 
Some examples are briefly discussed below (Gaussian, bi- 
modal and exponential distributions), to show the qual- 
itative differences between these computed pressure de- 
pendences and the experimentally observed features. 

For convenience, all computed values of 1/D e ff are dis- 
cussed in terms of P c ff = X> 3 / D e s (the pressure of freely 
diffusing 3 He gas required to yield a diffusion coefficient 
value equal to -D c ff)- Hereafter, the probing gas pres- 
sure is also simply designated by P 3 (equal to the actual 
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pressure for pure 3 He gas, or to the equivalent pressure 
P| q for a gas mixture as described in the text for the 
3 He — N 2 mixtures used in the experiments). 



1. Gaussian mfp distributions 

Figs. 11 and 12 display the results obtained for Gaus- 
sian mfp distributions characterized by mean values A 
and widths SX. 

A narrow Gaussian distribution leads to a variation of 
P e ff that is very similar to that described in Section V- 
A of the related article for a uniform medium with a 
single characteristic mfp A aoro equal to A . For instance, 
selected data from computations performed for various 
mean values Ao and constant width SX are presented in 
Fig. 11. The linear plot in Fig. 11 shows that, for each 
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FIG. 11: Broken lines: Pressure dependence of P e g =T>z/D e s 
[Eq. (S-2), with /(A) oc exp[-(A - \ {) ) 2 /SX 2 ]], computed for 
various mean values Ao (see legend) and <5A = 0.01 pm. Solid 
line: free gas data (P e s = f^)- 

Ao, P e ff actually scales linearly with P3 with an offset 
value that increases as Ao decreases. With 5X = 0.01 /zm, 
over the probed pressure range, all curves can be fit by 
straight lines with slopes equal to 1 and offset values 
Peff' (0) = 32?3 /vXq with an accuracy better than 50 ppm. 
The log plot in Fig. 11 therefore exhibits, for each value 
of A , a crossover from the free gas behavior (solid line) 
at high pressure, to a plateau at low pressure (where the 
mfp for free diffusion becomes so large that the effec- 
tive mfp is just that imposed by the aerogel, A ). This 
crossover occurs when the mfp A3 Ho for gas collisions be- 
comes comparable to Aq (i.e., around P3 = 32?3/i>Ao). 



The crossover region typically extends over one decade 
in pressure. 

The influence of the width <5A is depicted in Fig. 12, 
where Ao is kept constant and <5A is varied (Ao = 0.3 /zm; 
SX = 0.01 to 30 11m). The impact of the Gaussian width 
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FIG. 12: Broken lines: Pressure dependence of P c g — T>3/D c s 
for Gaussian mfp distributions in aerogel [Eq. (S-2), with 
/(A) a cxp[— (A — Ao) 2 /5A 2 ]], computed for various widths 
SX (see legend) and constant mean value A = 0.3 fim. Solid 
line: free gas diffusion (P c b = P3). Inset: Illustration of two 
extreme cases: broad (dashed line) and narrow (solid line) 
Gaussian mfp distributions. For the broad one, truncation 
at I = leads to an average mfp that is significantly higher 
(hence, to less restricted atomic diffusion). 

is quite small as long as SX ^ A , and corresponds to a 
slight change of the effective slope (linear plot in Fig. 12: 
dash, dot, and dash-dot lines). A significant change oc- 
curs when SX becomes comparable to Ao, due to trunca- 
tion at small SX. This situation is qualitatively depicted 
by the inset in Fig. 12 that displays Gaussian distribution 
profiles for SX <C A (solid line) and for SX > X (dashed 
line). For a further increase of the Gaussian width, the 
impact becomes more significant since much larger mfp's 
are reached (looser spatial confinement of the atoms), 
and a strong offset decrease is observed for SX 3> Ao 
(linear plot in Fig. 12: dash-dot-dot, short dash, and 
short dot lines). For instance, in spite of the short mfp 
mean value selected here (Ao = 0.3 /im), the free-diffusion 
limit would nearly be recovered over the experimentally 
probed pressure range (i.e., P3 > 10 mbars, and hence 
Agas < 40 /jm) for SX on the order of 30 /um (log plot in 
Fig. 12). 
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As a conclusion, whatever the choice of numerical pa- 
rameters, a Gaussian distribution fails to provide a suit- 
able pressure dependence to account for the experimental 
measurements. 



2. Bimodal mfp distributions 

Figs. 13 and 14 illustrate the expected behaviors for bi- 
modal distributions that are sums of two Gaussian func- 
tions with identical amplitudes. 

In Fig. 13, a narrow Gaussian distribution /W (Ao = 
50 /im, 5X = 0.1 /im) is combined with another narrow 
Gaussian distribution (mean: 10 /im, width: 0.1 /im). As 
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FIG. 13: Pressure dependence of P e ff obtained for a bimodal 
mfp distribution (solid line) that combines two narrow Gaus- 
sian distributions : Ao = 50 /im, <5A = 0.1 /jm (distribution 
dashed line), and A = 10 /im, SX = 0.1 /jm (dotted line), 
in log plot. Inset: linear plot, for the same distributions. 

expected, the dominant contribution arises from the dis- 
tribution that provides access to the largest mfp's. The 
relative difference between the bimodal distribution and 
the Gaussian distribution is noticeable at low pres- 
sure. The combination with the other Gaussian results 
in a significant offset increase (70% change), and a slight 
curve distortion below P3 = 100 mbars (best seen in the 
inset). The corresponding relative change in P c g values is 
plotted against gas pressure in Fig. 14 (solid line) . It be- 
comes negligible as soon as the free gas mfp gets smaller 
than the lowest bound of the mfp distribution in aerogel. 

Fig. 14 also displays results obtained for different val- 
ues of the width SX of the /W Gaussian distribution, 
that remains centered around A = 50 fim.. The relative 
impacts of the combination with the other Gaussian dis- 
tribution can be checked to become smaller as the width 
SX is increased, due to the reduced relative weight of this 
additional contribution (Fig. 14, broken lines). 

Whatever the combination of Gaussian distributions 
performed, bimodal distributions can absolutely not be 
adjusted to match the power-law pressure dependence 
that is experimentally observed. 
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FIG. 14: Pressure dependence of the relative changes in P ef i 
for bimodal mfp distributions obtained by combination of 
Gaussian distributions centered around Ao = 50 /im with 
various widths 5\ (providing the reference values P^g), and 
a fixed Gaussian distribution (mean: 10 /im, width: 0.1 /im). 
Solid line: SX = 0.1 /mi (value identical to that used in 
Fig. 13). Dashed line: SX — 1 /tm. Dotted line: SX = 10 /im. 



3. Exponential mfp distributions 

Fig. 15 displays the pressure dependences typically 
obtained for exponential distributions with characteris- 
tic mfp's A c ranging from 0.05 to 1 /im (identified by 
distinct line styles in the graph). P c g remains hnite at 



10000- 



-Q 

E 



5000- 



1000 



2000 



3000 



P„ (mbar) 



FIG. 15: Pressure dependence of P c s for exponential mfp 
distributions [Eq. (S-2), with /(A) oc exp( — A/A c )] for various 
characteristic mfp's A c : 0.05 /im (dashed line), 0.1 /tm (dotted 
line), 1 /tm (dash-dot line), and 10 /tm (dash-dot-dot line). 

very low pressures and approaches P c ff (0) = 3£>3/vA c , 
the expected value for an unbounded distribution. P c g 
also increases with pressure, but not linearly. The deriva- 
tive dP e g/dP3 smoothly varies between two asymptotic 
values: 2, at null pressure, and 1, at infinite pressure. 
For each P e g curve, the change of local slope essentially 
occurs over the pressure range where the mfp for free dif- 
fusion is larger than (or equal to) the characteristic scale 
A c introduced by the aerogel. Beyond this limit, the gas 
motion is hardly limited by the solid matrix. This nearly- 
free diffusion regime (linear pressure dependence with 
slope 1) is reached for P3 3> 40 mbars at A c = 10 /im, for 
instance. At A c = 0.05 /im, the probed pressure range 
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must be extended to P3 » 8 bars to observe it (data not 
shown in Fig. 15). 

When a lower bound I is arbitrarily set to the expo- 
nential mfp distribution, the results for P e g are strongly 
altered (Fig. 16). Since all mfp's smaller than I are re- 
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FIG. 16: Pressure dependence of P e e for truncated exponen- 
tial mfp distributions. Computations are performed for three 
values of the characteristic mfp A c : 0.05 /jm (dashed lines), 
0.2 /im (dotted lines), and 5 /jm (solid lines). For each A c 
value, three different lower bounds I are set to the mfp distri- 
bution: 0.001 /im (squares), 0.01 /im (circles), and 0.03 /im 
(triangles). 

moved from the distribution, a global reduction of effec- 
tive pressure may be expected (lesser gas confinement). 
Quantitatively, the new value expected at zero pressure 
is P e ff(0) = 3V 3 /v(X c + 1), and the relative impact of the 
truncation rapidly increases as l/X c increases. The trun- 
cation also induces a change of the P c g increase rate over 
the whole pressure range, related to the modification of 
the probability distribution. This effect depends on the 
lower bound values I, as shown in Fig. 16 for three expo- 
nential distributions (here also, line styles vary with char- 
acteristic mfp values A c ). For instance, detailed compari- 
son of the results obtained with lower bounds I = 0.01 /im 
(circles) and I = 0.03 /im (triangles), with reference to 
those obtained with I = 0.001 /im (squares) in Fig. 16 
yields the following results: 1/ The relative P e ff changes 
computed at fixed pressure, for P3 = 4.8 bars, amount to 
30%, 10% and 0.5% for A c = 0.05, 0.2, and 5 /im (respec- 
tively) with lower bound I = 0.01 /im. They drop to 14%, 
4% and 0.2% for A c = 0.05, 0.2, and 5 /im (respectively) 
with lower bound I = 0.03 /im. 2/ These P e ff changes de- 
crease slowly, but steadily, as the gas pressure increases. 
From Pj = 1 mbar to P 3 = 10 bars, the exhibited varia- 
tion is equal to 27% with lower bound I = 0.01 /mi, and 
drops to 14% with lower bound I = 0.03 /an, for the three 
A c values. 

The distribution may also be truncated at large mfp's 
to probe the sensitivity of the diffusion coefficient to 
such a limit. The influence of the upper bound L re- 



mains extremely limited due to the exponentially van- 
ishing statistical weight of the large aerogel mfp's, as 
long as the selected L value exceeds that of the char- 
acteristic mfp A c . Conversely, it is important for upper 
bounds L comparable to A c . The effective pressure at 
zero gas pressure depends, quite naturally, on the ex- 
ponential factor r = cxp(— L/A c ), and it is expected to 
vary as P cff (0) = (3V 3 /v\ c ) x [1 + rln(r) /(l - r)]" 1 . 
But P e ff changes must noticeably vary with gas pres- 
sure, being much larger at low pressures (restriction of 
atomic diffusive motion over large length scales has been 
lifted) than at high pressures (atom-atom scattering al- 
ready sets a small mfp limit). In Fig. 17, dotted lines 
show, for instance, the results obtained for truncated ex- 
ponential distributions with A c = 0.5 /im, and compara- 
ble upper mfp bounds L. Symbols are used to specify 
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FIG. 17: Pressure dependence of P e s for an exponential mfp 
distribution, with (broken lines) or without (solid line) trun- 
cation at large mfp's. The characteristic mfp of the distri- 
bution is A c = 0.2 /im. Data are computed for three up- 
per bound L values: 0.1 /jm (squares), 0.2 /im (circles), and 
0.5 /im (triangles). 

the selected L values, that vary from 0.1 to 0.5 /im (see 
caption). The P 3 = 1 mbar data differ from the zero 
pressure expected values P c ff(0) by less than 0.2%, for 
the three L bounds. With reference to the unbounded 
distribution data (Fig. 17, solid line), the P e ff increases 
due to the truncation vary from 4.35 (respectively: 2.39, 
and 1.29) at P3 = 1 mbar, to 1.48 (respectively: 1.19, 
and 1.04) at P3 = 10 bars, for L = 0.1 /im (respectively: 
L = 0.2 /im, and L — 0.5 /im). 

With or without truncation, exponential mfp distribu- 
tions fail to provide effective diffusion coefficients that 
behave like the experimental ones. 

In summary, for the three types of mfp distributions 
presented here, the computed pressure dependences ex- 
hibit features that do not correspond to the observed 
ones. A satisfactory description of the experimental dif- 
fusion data can only be obtained with the phcnomenolog- 
ical power-law distributions introduced and characterized 
in the related article (Section V-B). 
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